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ABSTRACT 



Aims. Chemical abundances in solar-type stars are still a much debated topic in many respects. In particular, planet-hosting stars are 
known to be metal-rich, but whether or not this peculiarity applies also to the chemical composition of the outer stellar atmospheres 
is still to be clarified. More in general, coronal and photospheric abundances in late-type stars appear to be different in many cases, 
but understanding how chemical stratification effects work in stellar atmospheres requires an observational base larger than currently 
available. 

Methods. We obtained XMM-Newton high-resolution X-ray spectra of t Bootis, a well known nearby star with a Jovian-mass close-in 
planet. We analyzed these data with the aim to perform a detailed line-based emission measure analysis and derive the abundances 
of individual elements in the corona with two different methods applied independently. We have compared the coronal abundances of 
T Bootis with published photospheric abundances based on high-resolution optical spectra and with those of other late-type stars with 
different magnetic activity levels, including the Sun. 

Results. We find that the two methods provide consistent results within the statistical uncertainties for both the emission measure 
distribution of the hot plasma and for the coronal abundances, with discrepancies at the Icr level limited to the amount of plasma 
at temperatures of 3-4 MK and to the O and Ni abundances. In both cases, the elements for which both coronal and photospheric 
measurements are available (C, N, O, Si, Fe, and Ni) result systematically less abundant in the corona by a factor 3 or more, with the 
exception of the coronal Ni abundance which is similar to the photospheric value. Comparison with other late-type stars of similar 
activity level shows that these coronal/photospheric abundance ratios are peculiar to t Bootis and possibly related to the characteristic 
over-metallicity of this planet-hosting star. 

Key words, stars: late-type - stars: atmospheres - stars: coronae - stars: individual: t Bootis - stars: abundances - X-rays: stars 



1. Introduction 

Chemical abundances in solar-type stars are still a much de- 
bated topic. Solar abundances have been derived both from 
optical/UV/X-ray spectroscopy and from c omposition measure- 
ments of solar wind particles or meteorites dFeldman & Laming, 
2000), but the most recent and sophisticated results are in stark 
conflict with models o f the solar interior tuned with helioseis- 
mology data ( Asplund et al.L 12009) . A second problematic point 
is the finding that photospheres of stars hosting close-in giant 
planets are known to be metal-rich, which can possibly be ex- 
plained by a higher probability of planet formation in high- 
metallicity birth clouds ( Gonzalezl I2006h . 

Another conundrum is offered by the differences found in the 
Sun and other solar-like stars between photospheric and coro- 
nal abundances. In this respect, stellar X-ray spectroscopy is a 
fundamental tool for the chemical analysis of stellar outer atmo- 
spheres, and the only method to determine the abundances of 
noble gases, like argon and neon, with no optical lines in photo- 
spheric spectra. For the solar corona, and in particular in long- 
lived coronal structures, the composition of the plasma appears 
enriched with elements having low first ionization potentials 
(FTP < 10 eV) by abo ut a factor 4, on average, with respect to 
photospheric values (iFeldman & Lamingl l2000l) . In other stars 
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a mo re complex behavior has been observed ("Gudel & Naz^ 
with a tendency for the low-FIP elements (including iron) 
to become depleted with respect to the high-FIP elements (neon 
in particular) in extremely active RS CVn-type and Algol-type 
binaries. Several theoretical explanations have been proposed, 
but our understanding is still largely driven by the observations, 
which indicate a dependence of the coronal/pho tospheric abun- 
dance ratios on the stellar ma gnetic activity level (iRobrade et al.L 
120081) or on the spectral type dWood & Linskvil2010l) ! 

The aim of the present work is to investigate the coro- 
nal/photospheric abundance problem in a metal-rich planet- 
hosting star To this aim we made use of detailed analyses of the 
photospheric composition, available in the literature, for stars of 
this class, and of X-ray data acquired on purpose with the XMM- 
Newton satellite. The comparison of our target with other solar- 
type stars with similar levels of magnetic activity provides new 
information on the relative importance of the photospheric com- 
position in the abundance stratification mechanism(s) operating 
in stellar atmospheres. In Sect.|2]we present our sample star, and 
in Sect.|3]its photospheric chemical composition. The analysis of 
the X-ray spectra, including the derivation of the plasma emis- 
sion measure distribution (HMD) vs. temperature and coronal 
abundances is described in Sect. |4] Sections|5]and|6]are devoted 
to discussion and conclusions. 
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2. Why rBootis 

rBootis (HR 5185, HD 120136) was selected as part of a project 
aimed at studying the coronal abundances of late-type dwarfs 
classified as super-metal-rich stars, and sufficiently X-ray bright 
to obtain good X-ray spectra with the XMM-Newton Reflection 
Grating Spectrometers (RGSl and RGS2). 

r Boo is a 1 .3 M© F7 V star (i Perrin et allll977l:lBaines et all 

I2008h . located at 15.6pc from the Sun. It harbors a planet with 
a mass msin/ k 3.9 Mj, discovered bv Butler et alj dT997 ^ in a 
3.31 day period. The stellar age was estimated by several authors 
with different techniques (in particular, position in the H-R dia- 
gram and chromospheric emission level) and constrained within 
the range 1-3 Gyr (Takeda et al., 200 ^. How ever, based on its 
X-ray luminosity (Lx ~ lO^'^erg s"', Schmitt et al., 19 85|), the 
age o f r Boo could be as young as 0.4 Gyr (.Sanz-Forcada et al.L 
120101) . 

The rotatio n characterist i cs of r Boo were determ i ned b y 
iReinersI (l2006l) . ICatalaetalJ (l2007h . and iDonati et al.l (l2008h . 
The photospheric line profiles are reasonably well fitted as- 
suming an homogeneous surface rotational velocity vsin/ - 
15 + 1 kms"' and a turbulent velocity of 5.5 km s"' . Fourier anal- 
ysis of the same data was employed to establish a surface differ- 
ential rotation of ~ 20%, which implies a rotation period at the 
equator of 3.0 d, and 3.7-3.9 d at the pole, assuming a rotation 
axis inclined at an angle / = 40° with respect to the line of sight. 
Hence, rBoo is considered the first convincing case of a star 
whose rotation (at intermediate latitudes) is synchronized with 
the orbital motion of its close-in giant planet. 

ICatala et al.1 (l2007h and lDonati etal] (l2008l) also studied the 
strength, topology, and time change of the surface magnetic field 
of r Boo. At large spatial scales, the field comprises a dominant 
poloidal component, more complex than a dipole, and a small 
toroidal component. Field strengths of up to 10 G were found, 
and the overall polarity turned out to be reversed between two 
observations taken about one yea r apart. Mos t recen tly, a second 
polarity reversal was reported bv lFares et al ] (l2009h . suggestmg 
a magnetic cycle of about two years. 

Finally, rBoo is also one of the stars recently investi- 
gated to search for e vide nce of star-p l anet i nteraction effects. 
IShkolnik et all (l2008h and lWalker et al.l (1200 8!) presented studies 
of the photospheric (white light) and chromospheric (Ca ii K line 
core) variability of r Boo over a time range of a few years, sug- 
gesting the presence of a persistent dark spot and enhanced chro- 
mospheric emission synchronized with the planet orbital period, 
but occurring at phase ~ 0.8 (i.e. preceding the sub-planetary 
stellar longitude, = 0, by ~ 70°). These results are suggestive 
of variability caused by an active region induced by a star-planet 
magnetic interaction, but the signal appears marginally signifi- 
cant, an d further studies ar e required to confirm this conclusion 
(see also lFares et alil2009l) . 

In summary, rBoo is a relatively massive star with a shal- 
low convection zone and a fairly massive planet in a tight orbit. 
The magnetic activity of r Boo is moderate and typical of mid-F 
dwarfs, with little if any modulation at either the orbital or rota- 
tion period. The near synchronization between the star's rotation 
and the planetary orbit might imply little or no stress of the mag- 
netic field lines that connect the two components of the systems, 
and hence a low planet-induced activity enhancement. In this re- 
spect, r Boo is an ideal target for studying element abundances 
in the quiescent corona of a star with a photospheric metallicity 
significantly exceeding the solar one. 



Table 1. Adopted photospheric abundances of Tau Boo 



Element [Ab/H] 





0.329 


± 0.084 


[N/H] 


0.2 


±0.1 


[O/H] 


0.300 


± 0.055 


[Na/H] 


0.463 


± 0.096 


[Al/H] 


0.398 


±0.100 


[Si/H] 


0.333 


± 0.073 


[Ca/H] 


0.161 


±0.104 


[Fe/H] 


0.332 


± 0.079 


[Ni/H] 


0.217 


±0.149 



Notes. Al l abundances a r e from iGonzalez & Law j d2007t) . except 
[N/H] from lEcuvillon et all d2004ah . 

3. Photospheric abundances 

The composition of Tao Boo's photosphere was determined by 
several authors in the context of many studies carried out to 
quantify and explain the peculiar metallicities of stars with plan- 
ets (SWPs; see Gonzalez, 2006, for a review). 

We sel ected r Boo o riginally from the extensive com- 
pilation of ^ Taylor! (11994'), who has critically examined the 
high-resolution spectroscopic abundance determinations of F, 
G, and K stars available in literature. rBoo was noted as 
a metal-rich star according to the criterion [Fe/H] > 0.2. 
Subsequently, the high metallicity of our target was confirmed 
by several accurate abundance analyses, the most recent per- 
formed b y Ecuvillon et al. (2004b), Valenti & Fischer (2005L 
Takeda & Honda (2005L iGilU et all (120061) . iLuck & HeitS 
(2006) . and Takeda (200% . 

Most of the abo v e stud ies were considered in the work 
by iGonzalez & Lawsl (l2007h . who have presented the chemical 
abundances for 18 elements in 31 SWPs and compared them 
with a sample of nearby stars without detected planets. In order 
to cope with differences in the zero-point, effective temperature 
and abundance scales among the various studies, these authors 
have applied statistical corrections to abundance measurements 
for individual ele ments reported by each author For these cor- 
rections, they set lGilli et alj (l2006l) as t he referenc J| fo r all the 
ele ments, except fo r C and O, for which lEcuvillon et al.l ( l2004bl) 
and iLuck & Heitej (l2006h . respectively, were adopted. Finally, 
they computed average abundance values for each element, and 
standard deviations from individual uncertainties summed in 
quadrature. Note that this procedure yields conservative (slightly 
overestimated) uncertainties on the final values, larger than the 
scatter between individual measurements. The abundance values 
of rBoo computed by Gonz alez & La ws (2007) are reported in 
Table [1] wh ere we also added the si ngle nitrogen abundance de- 
termined bv lEcuvillon et alj(l2004ah . These are the photospheric 
abundances we adopted as reference for comparison with the 
coronal abundances derived by us. 



In turn, the stellar abundances in iGilli et alJ 12003) are give with 
respec t to the solar photospheric abundances bv lAnders & Grevessd 
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Table 2. Line identifications and fluxes in the Tau Boo spectrum. 
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Notes. Observed minus predicted line flux, divided by flux error. Flag= 1 if line was employed for EMD reconstruction with Method 1 only, 
Flag=2 if Method 2 only, Flag=3 if employed by both methods. 
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Fig. 1. X-ray light curve of rboo obtained with EPIC (summed 
PN and MOS broad-band photometry). Time bins are 600 s long. 
The upper X-axis scales indicate the planetary orbital phase at 
the time of the observation. 

4. X-ray emission 

tBoo was observed on June 24, 2003, for a duration of about 
65 ks, with the thick filter in front of the EPIC detectors, the 
MOS in LARGE WINDOW mode, and the PN in full frame mode. 
After screening for high-background intervals, the useful expo- 
sure time was reduced to about 50 ks, during which the source 
yielded a mean count rate of 0.3 s ' in the MOS, and 1.3 s ' in 
the PN (0.3-5keV band). 

Figure [T] shows the X-ray light curve of t Boo with a bin 
size of 600 s, just characterized by some low-level variability, 
typical of solar-t ype stars with mod erate activity. According to 
the ephemeris bv lCatala et"an (l2007h 

r = HJD 2453450.984-1- 3.31245 £ (1) 

during the 16 hours of the XMM-Newton observation the planet 
orbital phase was between 0.89 and 1.13, i.e. the planet was near 
the oppositiorQ (occulted by the star). 

The X-ray spectral analysis was carried out in several steps: 
first, a preliminary fit to the EPIC data with three-component 
thermal models was performed to determine the level of contin- 
uum emission and to guess the plasma metallicity; these results 
were employed to proceed with the measurement of emission 
lines in the RGS high-resolution spectra, and eventually for a 
detailed line-based emission measure analysis, including an es- 
timate of chemical abundances for a number of elements with 
measurable spectral signatures (in particular, C, N, O, Ne, Fe, 
and Ni); the plasma emission measure distribution and avail- 
able element abundances were employed subsequently to de- 
termine high-temperature components and other element abun- 
dances (Mg, Si, S, Ar) by again fitting the EPIC spectra. Below 
we explain these steps more in detail. 

4.1. RGS data analysis 

The background-subtracted RGS 1 and RGS2 spectra of r Boo 
contain ~ 4000 and ~ 5100 net counts, respectively (Fig.|2|i. 



^ Note that the zero-point phase in Catala et al. (2007) is shifted by 
half a period w ith respect to the ephemeris in Wa lker et al.. (2008.) and in 



The analysis of X-ray emission line spectra is notoriously 
difficult, and the measurement of element abundances can be af- 
fected by systematic uncertainties that are difficult to quantify 
because of the assumed plasma emissivity model, the identi- 
fication and measurement of relevant lines, and the algorithm 
adopted for the emission measure analysis (Maggio, 2007). For 
this reason, we employed two different approaches for the anal- 
ysis of the XMM-Newton grating spectra of r Boo, carried out 
in a completely independent manner. 



4.1.1. Method 1 



IShkolnik et al.l (2008) (based on older photometric data), where phase 
is with the planet in front of the star 



The fi rst approach is the one discussed in detail in 'Scel si et aP 
(|2004|) together with a study of its accuracy; here we limit our- 
selves to report the main points of this iterative method. 

We employed the s oftware package PINTofALE V2.0 
dKashvap & Drake*, '2000") and, in part, also XSPEC, and used 
the As trop hysical Plasma Emission Database ( APE D/ATOMDB 
y 1.3, [Smith etali 12001). which includes the iMazzotta et all 
( Il998h ionization equilibrium. 

Background-subtracted RGSl and RGS2 spectra were first 
re-binned by a factor of 2 and then co-added for the identifi- 
cation of the strongest emission lines and the measurement of 
their fluxes. In the rebinned RGS1+RGS2 spectrum we iden- 
tified about 40 lines from Fe xvii-xxi, Ne ix-x, O vii-viii, N vii, 
C VI, and Ni xix-xx ions. After this step the method proceeds in 
an iterative way. 

To measure line fluxes we adopted a Lorentzian line pro- 
file and initially assumed the continuum level evaluated from 
a 3-T model best-fitting the pn spectrum. Indeed, the wide line 
wings make it impossible to determine the true source contin- 
uum below ~ 17 A directly from the RGS data, especially in the 
~ 10 - 17 A range where the spectrum is dominated by many 
strong overlapping lines. The Mgxi triplet was measured as a 
single feature by counting the photons within the 9.10 - 9.40 A 
range in the RGS spectrum rebinned by a factor of 5, and sub- 
tracting the assumed continuum contribution in the same wave- 
length region. We also fixed A^h = 3 x 10'^ cm"^, compatible 
with the short distance of t Boo (15.6 pc). The results are listed 
in Table |2l 

To reconstruct the emission measure distribution (EMD, de- 
fined as EM(T) = J^^ riHtii-dV [cm ]) and to estimate element 
abundances, we selected a set of 19 lines among those identi- 
fied that had reliable flux measurements and theoretical emissiv- 
ities. Most of them are blended with other lines, so the measured 
line flux is actually the sum of the contributions of a number 
of atomic transitions. Accordingly, the "effective emissivity" of 
each line blend must be evaluated as the sum of the emissivities 
of the lines that mostly contribute to that spectral feature. If the 
blend includes lines from different elements, the effective emis- 
sivity needs to be weighted by their relative abundances. Since 
these abundances are not available at the beginning of the anal- 
ysis, we initially selected only iron lines not blended with lines 
of other elements, and the latter are included one at a time in the 
following stages. 

We performed the EMD reconstruction with the Markov- 
Chain Monte Carlo (MCMC) method by iKashvap & Drake! 
(1998). This method yields a volume emission measure distri- 
bution, EM{Tk) - dem(Tk) AlogT, and related statistical un- 
certainties AEM(Tk), where demij) - n^dV/dlogr is the 
differential emission measure of an optically thin plasma and 
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Fig. 2. Summed RGSl and RGS2 spectra of rBootis (histogram) with the synthetic model spectrum derived from the analysis 
results obtained with Method 1. Labels identify some of the strongest lines. 



A log r = 0.1 is a constant bin siz^B The method also provides 
estimates of element abundances scaled by the iron abundance, 
with their statistical uncertainties. The iron abundance is esti- 
mated by computing the emission measure distribution assum- 
ing different metallicities and by comparing the synthetic con- 
tinuum with the observed one at /I > 20 A in the RGS spectrum, 
which is a spectral region free of strong overlapping emission 
lines. Finally, we checked the solution obtained with the MCMC 
by comparing (i) the line fluxes predicted from our solution with 
the measured ones and (ii) the whole model spectrum with the 
observed one (Fig.|2]l. 

The procedure described above is repeated two or three times 
to check line identifications and to ensure consistency between 
the continuum level assumed for flux measurements and the con- 
tinuum predicted by the EMD. Indeed, the continuum is adjusted 
at each iteration because it may become different from that pre- 
dicted by the 3 - T model best-fitting the pn spectrum, which is 
adopted as the initial guess. This procedure ensures that possible 
cross-calibration offsets between pn and RGS do not affect the 
final EMD. 

The final EM(T) values are reported in Table|3]with the coro- 
nal abundances of individual elements. The EMD is also shown 
in Fig. [3] with the distribution obtained with Method 2. 
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Fig. 3. Emission measure distribution of t Boo based on the 
XMM-Newton RGS sp ectra. The th i ck (re d) histogram was ob- 
tained with Method 1 (IScelsi et all l2004. the th in (blue) his- 
togram with Method 2 (ISanz-Forcada et al.Ll2003h . In flie latter, 
dot-fined bins are more uncertain (no associated error bar). 



4.1.2. Method 2 

The second approach for constructing the E MD vs. temperature 
follow s the line -based analysis described in lSanz-Forcada et al.l 
(12003 ) and references therein. In short, individual line fluxes are 
measured and those with a signal-to-noise ratio S /N > 3^ are 
then compared with the theoretical values obtained from a trial 
EMD, combined with the same atomic emission model used for 



^ We also performed the EMD analysis with a coarser A log T = 0.2 
temperature grid, but obtained consistent results within statistical un- 
certainties. 



Method 1 . The comparison of measured and model line fluxes 
results in an improved EMD that is used again to produce new 
model line fluxes. 

In Method 2 we employed 27 lines. The first iteration is 
based only on Fe lines, either blended or unblended, assuming 
solar abundances. The successive iterations are used to calculate 
the abundances of the different elements, that are required to re- 
fine in particular the line predictions for all the blends. Because 
the lines corresponding to the various elements are formed in dif- 
ferent temperature ranges, we take the caution to progressively 
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Table 3. Volume emission measure distribution and coronal 
abundances of Tau Boo 
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Notes. All abundances are relative to the solar photospheric val- 
ues bv lAnder~& Grevesse. L1989.) . Values with no error were fixed (see 
text). 



add lines of elements with formation temperatures that overlap 
the EMD computed at each iteration. 

The method converges in a few iterations, and is therefore 
self-consistent. This process yields an EMD solution that is not 
unique, but reliably matches the model fluxes to the observed 
ones, and it presumably describes the actual amount of plasma 
at diff'erent temperatures in corona. 

Finally, error bars are calculated using a Monte Carlo method 
that seeks the best solution for d ifferent line fluxes within their 
1-cr errors (see lSanz-ForcadaTT al.. 2003, for more details). 

In summary, the two methods differ in the choice of the line 
set employed for the EMD reconstruction, in the algorithm to 
search for the best EMD solution, and in the evaluation of its 
confidence range. 

4.2. EPIC spectral analysis 

Typical EPIC spectra can usually be fitted with thermal mod- 
els with few (2 - 4) isothermal components. This is an approx- 
imation to the actual emission measure distribution of the coro- 
nal plasma, which provides a reasonably good quality of the fit 
to individual CCD-resolution X-ray spectra (MOS or PN); in- 
stead, systematic errors and cross-calibration uncertainties make 
the joint analysis of all the available spectra (EPIC and RGS) 
cumbersome, and yield best-fit values formally too high to be 
statistically acceptable. This is also the case for t Boo. 

The simplest multi-component thermal model that provides 
an acceptable fit to the EPIC spectra is a three-temperature 
model with plasma at 2MK, 4MK, and 8MK and a volume 
emission measure in the proportions 1:6:1 (see Table [3] and Fig. 



T Boo, PN 




0.5 1 2 5 10 

Energy [keV] 



Fig. 4. Upper panel: PN observed spectrum (black histogram 
with error bars)) and best-fit model (in red). Lower panel: resid- 
uals in sigma units. 



|4]i. This model fitting confirms the low interstellar absorption 
assumed in the analysis of the RGS spectra (Sect. 14.1.11 ). and 
yields an independent measure of the abundances of O, Ne, 
Mg, Si, and Fe. Note that Mg and Si abundances can be ef- 
fectively constrained by this fitting, because the EPIC spectra 
clearly show line complexes caused by Mgxi (1.33-I.35keV), 
Mgxii (1.47keV), and Sixiii (1.84-1.86keV). Abundances of 
other elements cannot be constrained by the EPIC spectra, and 
they have been fixed either to the values derived from the analy- 
sis of the RGS spectra, or to the iron abundance; in this way, the 
best-fit model also reproduces sufficiently well the weak spectral 
feature visible at the expected position of the S xv triplet (2.43- 
2.46 keV). 

5. Discussion 

5. 1 . Abundance depletion in the corona of t Boo 

The central result of our analysis is illustrated in Fig.|5j the coro- 
nal abundances of our target are all systematically lower than the 
photospheric ones, with the only possible exception of the Ni. In 
this respect, there is agreement on the abundances derived with 
any of the methods we employed for all elements but oxygen 
and nickel, within the statistical uncertainties. The 2cr difference 
in the former case can be easily explained by inspection of Fig. 
[3] where the emission measure derived with Method 1 in the 
temperature range log T = 6.3-6.5 is systematically higher than 
indicated by Method 2: since much of the Ovii and Oviii line 
fluxes are produced in this range, an uncertainty in the emission 
measure implies a systematic error in the abundance of this el- 
ement. The case of the Ni abundance is less obvious, but likely 
because in APED/ATOMDB V1.3 the Ni atomic model is less 
accurate than for other elements, and the Ni lines are in general 
relatively faint and hence affected by larger measurement uncer- 
tainties. 

Regardless of the method, our results show that the coronal 
iron abundance is about half the solar photospheric value, while 
the stellar photosphere is twice the same value. The statistical 
uncertainty on the ratio of the photospheric to coronal iron abun- 
dance is very small: 0.24 + 0.09 with Method 1, and 0.23 + 0.06 
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Fig. 5. Photospheric and coronal abundances of r Boo, sorted by 
increasing Fi rst Ionization Potent i al. Th e reference solar abun- 
dances are bv lAnders & Grevessd (Il989h . 



with Method 2 (uncertainties computed in quadrature), but an 
additional (systematic) source of uncertainty on the absolute 
coronal abundances derived from high-resolution X-ray spec- 
tra can be caused by the determination of the continuum emis- 
sion level. However, the completeness and accuracy of atomic 
databases for the synthesis of multi-temperature model spec- 
tra w as significantly improved in the last decade (iBrickhousel 
I2008h . and the results of our RGS data analysis are confirmed by 
the spectral fitting of the EPIC data as well. 

We conclude that our results are sufficiently robust to state 
that the X-ray emitting plasma of r Boo is significantly depleted 
in metals compared with the photosphere. Although the deter- 
mination of the absolute iron abundance in the corona may be 
affected by an uncertainty larger than indicated by our analysis, 
the pattern of abundances scaled by Fe is well constrained and 
the available RGS and EPIC spectra exclude a plasma metallic- 
ity matching the stellar photospheric composition. 

Moreover, we find little if any dependence of the coronal 
abundances on the first ionization potential (FIP), a result which 
holds true also considering the coronal vs. photospheric abun- 
dance ratios (Fig. |6]l, rather than values with respect to the so- 
lar chemical composition. In this respect, our results would not 
change after the po ssible choice of sola r element abundances 
more recent than the lAnders & Grevessd (fl989) reference scale, 
which we preferred for consistency with the majority of the pre- 
vious abundance studies. 



5.2. T Boo: a medium-activity coronal source 

In order to put the coronal abundances of t Boo in context. 
Fig. Q shows a comparison of its EMD distribution vs. temper- 
ature ( Method 1) with th ose of other fo ur G-K stars: the Su n 
(G2V, iPeres et al.Ll200(]h. Ceti (GSV , iTelleschi et all l2005h . 
6Eri (K2V , 'Sanz-ForcadaetaLll2QQl, and HD 283572 (G2IV, 
IScelsi et a l., 2005). 

The relative amount of plasma at different temperatures is a 
way to judge the level of stellar magnetic activity, and it carries 
more information than a simple comparison of broad-band X- 
ray luminosities. Among the stars we selected, the Sun is clearly 
the one with the lowest average temperature, with an EMD that 
peaks just below 2 MK. At the other extreme, HD 283572 is a 
young weak-lined T Tauri star in the Taurus-Auriga region with 
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Fig. 6. Upper: Coronal to photospheric abundance ratios for the 
Sun (black squares), t Boo (red circles and filled diamonds), and 
two other intermediate-activity stars (magenta empty diamonds 
for /<•' Cet and brown triangles for eEri). Elements are sorted as 
in Fig. |5] Lower: Same plot but with all values divided by the 
Fe abundance in a logarithmic scale. The FIP bias in the solar 
corona is indicated by the shaded (orange) strip, and for t Boo 
the abundances obtained with Method 2 are not shown for clar- 
ity. 



a much hotter and X-ray bright corona at 10-20MK. The other 
three stars were selected instead for the similarity of their EMD 
with that of r Boo. eEri is also known to host a Jovian-mass 
planet (m sin / 1 .6 Mj) at a distance of 3.4 AU, while /c' Cet is 
a putative single star. 

In passing we note that tBoo also has a s tellar compan- 
ion, GJ 527 B, with spectral type M2 (V=ll.l) jPatience et all 
2002). We evaluated the possible contamination to the observed 
X-ray flux by deriving the B-V color and the bolometric cor- 
rec tion from the s pectral type and adopting the transformations 
bv ' Flowed (fT996h . and thus obtaining for the bolometric lumi- 
nosity logLboi = 30.96 ergs"'. Assuming, conservatively, that 
the GJ 527 B X-ray emission is at the saturation level, we get 
Lx = lO^-'Lboi ^ 10^^ ergs"', which is about a factor 10 lower 
than the X-ray luminosity of tBoo (logLx = 28.94 ergs ', in 
the 0. l-2.4keV band). We conclude that the contamination from 
this unresolved companion is negligible. 

We took the photos pheric abundances for eEri from 
iGonzalez & Lawsl (l2007h . the same source as for rBoo, 
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Fig. 7. Emission measure distribution of t Boo (Method 1) com- 
pared with those of other G- or K-type stars with different ac- 
tivity levels, including the quiescent solar corona (see text for 
references). 

while for Cet we adopt ed the photospheric composition by 
lAUende Prieto et a n (12004 . which was determined in a strictly 
differential way with respect to the Sun. The coronal abundances 
of the sa me two stars were d erived by Sanz-Forcada et al. ( 2004) 
and by Telles chi et al.l ("2005 respectively, with an analysis of 
XMM-Newton data similar to ours. We then computed the coro- 
nal/photospheric abundance ratios of all the solar-type dwarfs in 
our sample, excluding HD 283572 because of its peculiar evolu- 
tionary phase and extreme activity level. 

The pattern of abundance ratios vs. FIP (Fig.|6]l appears very 
similar for tBoo and eEri, and significantly different from the 
solar case, which is charact erized by the well-known FIP effect 
('Feldman & Lamingl |2000'). Finally, a:' Cet displays an over- 
abundance of low-FIP elements (Mg and Fe) in the corona, al- 
though with relatively large uncertainties on the abundance ra- 
tios, but high-FIP elements (C and O) are underabundant with 
respect to the photospheric composition, which is not the case in 
the solar corona. In the right panel of Fig.|6]we also show a com- 
parison of all the coronal/photospheric abundance ratios scaled 
by the iron abundances, because relative abundances are con- 
sidered more robust than absolute abundances. This plot clearly 
shows that in tBoo and eEri there is a smaller (if any) depen- 
dence of the abundance ratios on FIP than in the solar case. 

The above comparison suggests that the outer atmospheres 
of some intermediate activity stars, like tBoo and eEri are 
not affected by a FlP-dependent chemical stratification as pro- 
nounced as in the solar case. A closer look at the left panel in 
Fig.|6]shows however that only for t Boo the coronal abundances 
are significantly and uniformly lower than the photospheric ones 
(with the possible exception of Ni). 

Common understanding is that the abundance patterns in 
low- and high-activity stars are different, with low-activity stars 
showing an overabundance of low-FIP elements in the corona 
(positive FIP effect), while high-activity stars are character- 
ized b y the reverse behavior (inverse FIP effect) dAudard et al.L 
l2003l) . Following this activity-dependence scenario (see e.g. 

For k"^ Cet we adopted the coronal abundances de rived with the 
line-b ased method an d the APEC plasma emissiv ities in iTelleschi et alJ 
( l2005h . scaled to the I Anders & Grevessd ( 119891) solar reference abun- 
dances. 



iGiidel & Nazi. l2009h . an abundance ratio Ne/Fe ~ 0.25 (with 
respect to the photospheric value) is expected for Sun-like stars, 
while Ne/Fe ~ 5-10 in very active stars. An intermediate activ- 
ity star with an EMD similar to those of tBoo, eEri, and /c' Cet, 
should display Ne/Fe ~ 1 . 

However, this scenario was constructed assuming stellar 
photospheric abundances similar to the solar ones - a com- 
mon practice for many years - which has recently been chal- 
lenged by studies where coronal abundances are compared to 
accurate measurements of stellar photospheric abundances: in 
most cases, this properly consistent comparison has shown little 
if any signific ant difference betwe en photospheric and coronal 
compositions (ISanz-Forcada et al.L 12004: Maggio et al., 20071 : 
ISanz-Forcada et al.Ll2009l )~ 

The novel result of our work is that tBoo does not show 
any evident FIP bias, but that on the other hand the coronal 
abundances are all significantly lower with respect to the pho- 
tospheric ones by factors of 3 to 9 (taking into account the un- 
certainties, and with the exception of Ni). Hence, the case of 
tBoo demonstrates that coronal and photospheric abundances 
may differ in a peculiar FlP-independent way, which was never 
reported this clearly before the present study. 

This result represents a new challenge for theoretical models 
that try to explain element fractionation in stellar outer atmo- 
spheres. In spite of many attempts in the past, the only model 
that is capable of producing eithe r ji pos i tive or an inv erse FIP 
effect is the one proposed by Lamingl (|2004 l2009h . In this 
model charged ions are subject to a ponderomotive force caused 
by Alfven waves propagating through the chromosphere, which 
may originate either in the corona or in the convection zone. 
Depending on whether the waves are reflected or not at the base 
of coronal magnetic loops, the ponderomotive force can point 
upward or downward, producing opposite effects. Low-FIP ele- 
ments are primarily affected because they are nearly fully ion- 
ized in the upper chromosphere, whereas high-FIP elements are 
partially neutral. In any case, the element fractionation is ex- 
pected to show a FlP-dependence, which is apparently lacking 
in the outer atmosphere of t Boo. 

6. Summary and conclusions 

We performed a detailed comparison of the photospheric and 
coronal abundances of the planet-hosting star r Boo, making use 
of the best optical and X-ray spectroscopic data available to date. 
We analyzed high-resolution spectra obtained with the reflec- 
tion gratings (RGS) on-board the XMM-Newton satellite, and 
also medium-resolution (EPIC) spectra with much higher pho- 
ton counting statistics. The emission measure distribution vs. 
temperature and element abundances of the coronal plasma were 
derived with two independent methods from the RGS data, and 
confirmed with the analysis of the EPIC spectra. 

In spite of some 2cr difference in the EMDs derived with 
the two methods, a robust and consistent result is that the coro- 
nal abundances of tBoo are significantly lower than the photo- 
spheric ones, by about a factor 4, on average. The largest uncer- 
tainty is for the oxygen abundance, whose underabundance fac- 
tor is in the range 4-9. We do not see any way to bring coronal 
and photospheric abundances into agreement, unless the coro- 
nal Fe abundance is 4 times larger than the present best estimate 
(s! 0.5 the solar one), but such a systematic error is very unlikely 
with data of the present quality and with the state-of-the-art anal- 
ysis methods we employed. 

Assuming solar abundances as reference, the pattern of coro- 
nal abundances vs. FIP is nearly flat, and this behavior is also 



8 



A. Maggio et al.: Photospheric and coronal abundances of rBootis 



confirmed when the actual tBoo photospheric abundances are 
considered. No important FlP-dependent element stratification 
mechanism appears to be at work in the outer atmosphere of 
T Boo, but the metallicity in the corona is very different from the 
photospheric one. It is unclear at the moment whether this is a 
peculiarity of metal-rich stars. 

From a complementary point of view, while the coronal char- 
acteristics of rBoo including absolute abundances are similar 
to other intermediate-activity stars, its photospheric abundances 
are recognized to be quite large for Galactic old-disk standards 
and linked to its being a star hosting a close-in giant planet. Our 
results challenge the present understanding of chemical strati- 
fication effects in the atmospheres of solar-type stars, and the 
possible role of stellar activity in the underlying physical mech- 
anism. 
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